A necessary and sufficient condition for optimal decompositions 
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An important measure of bipartite entanglement is the entanglement of formation, which is defined 
as the minimum average pure state entanglement of all decompositions realizing a given state. 
A decomposition which achieves this minimum is called an optimal decomposition. However, as 
for the entanglement of formation, there is not much known about the structure of such optimal 
decompositions, except for some special cases, like states of two qubits or isotropic states. Here we 
present a necessary and sufficient condition for a set of pure states of a finite dimensional bipartite 
systems to form an optimal decomposition. This condition is well suited to treat the question, 
whether the entanglement of formation is additive or not. 

PACS numbers: 03.67.-a 

One of the main ingredients in quantum information 
theory is entanglement. Therefore, there have been many 
attempts to understand this property better. While in 
the case of bi-partite pure states, entanglement is well 
understood and uniquely quantified under some general 
assumptions ||^ this is not the case for mixed states. 
There have been many different entanglement measures 
proposed for such states, one of which is the entanglement 
of formation |^, defined as 

Ef{p) min{^p,i?(|V'.»| ^kIV'.XV'zI = p}, (1) 

i i 

where E denotes the entanglement measure on pure 
states, 



tr (T log2 cr, a — trs 



Closely related to the entanglement of formation are so 
called optimal decompositions of a state, which are those 
decompositions {{pi, that achieve the minimum in 

eq. (|l|). However, for all cases except (g) and 
some highly symmetric states in higher dimensions |^ , it 
is not known, how to compute an optimal decomposition 
for a given state nor the entanglement of formation itself. 

One property of the entanglement of formation is be- 
ing afhne on the convex set of states, generated by con- 
vex linear combinations of pure states from an optimal 
decomposition In other words, if the decomposition 
{{pi, \ipi)}i'=i with Pi > is optimal, then the decompo- 
sition {{qi, where the Qi form an arbitrary prob- 
ability distribution is also optimal. The optimality of a 
decomposition thus only depends on the states and 
not on the corresponding probabilities. This is the jus- 
tification to talk about the optimality of a set of pure 
states by which we mean the optimality of the 
decompositions {(pi, |■^/'^))}^• 
A question concerning the entanglement of formation, 
which is widely believed to be true but not yet answered, 
is the question of additivity, i.e. whether Ef(p (g) a) ~ 
Ef{p) + Ef{a) holds true for all p and a. This question 
can be directly reformulated in terms of decompositions: 



Given the optimal decompositions {{pi, of p and 

iilj' of '^^ is tcnsorproduct decomposition 

{{piqj, which is clearly a decomposition 

of p(g) a, also an optimal one? In terms of optimal sets, 
this is equivalent to the question, whether the set of pure 
states {IV'i) ® is optimal given that 

and {\(l>j)}J^i are optimal. To address this question we 
present a necessary and sufficient condition, to decide 
whether a given set of pure states is optimal. 

Before we continue, we need to fix some notations. A 
decomposition of a state p is normally defined as a set of 
normalized vectors with appropriate probabilities pi 
such that p — X]i-P«IV'i)(V'j|- However, for computational 
and notational purposes it is often convenient to define 
a decomposition to be a set of non normalized vectors, 
whose squared norms correspond to the probabilities, i.e. 
we make the transition {{pi, \ipi))}i — > {I'/'i)}! with = 
^/pilV'i)- To distinguish non normalized from normalized 
vectors we mark them with a tilde. In this notation we 
have p = J2i \ 4>i){(i^i \ and 

Efip) = min {Y.£m)\J2 i^^x-^'i = p}^ 



where we have introduced the homogeneously extended 
pure state entanglement 



£m ■.= UfE{^) = -tva\og, 



tr a 



with a = ti's 

To find all possible decompositions {{qj, \4>j))}jLi of 
a state p from a given decomposition {(ft, ( 
e.g. the spectral decomposition) we will make repeated 
use of a theorem by Wootters et al. [||. This theo- 
rem states, that two decompositions {{qj, \4>j))}jLi and 
{{pi, generate the same state if and only if they 

are related by a right unitary matrix U, i.e. a matrix 
obeying WU = 1, in the following way: 



(2) 



2 



or in terms of non normalized decompositions 

n 



(3) 



We now state our main result, a necessary and suffi- 
cient condition for the optimality of a set of pure states 
on a bipartite finite dimensional Hilbert space: 

Theorem: Let \ipi) S Ha^^'Hb, i — 1,... ,n and 
(Jij :— trs |. The set of pure states {|V'i)}"=i is 

optimal if and only if 

n n 

c,\%l>i)) > - Re ( ^ QCj tr[cry loga an]) (4) 

i=l — 1 

holds true for arbitrary Ci € C 

Proof. First we prove sufficiency: Let {|V'i)}"=i be a set 
of pure states for which eq. (^) is satisfied. Choos- 
ing an arbitrary probability distribution {pijj'Li we can 
form a decomposition {(pi, representing the 

state p = Y^^=iPi\i^i){'^i\- Every other decomposition 
{('Zji l0j))}j^Li of p is obtained via a right unitary matrix 
f/by 



i=l 

From eq. we get 

n 

n 

> -Rc(^ X UjiUjk\/PiPk tr[o-ifc l0g2 CTr 



with Uife = tr Summing over j and using the 

right unitarity of U we arrive at 



> 



- Rc^Pi tr[CTii log2 an] = ^pii^dV'i)) 



which proves the optimality of the decomposition 
{(Pi ^^"i therefore the optimafity of the set 

Next we prove necessity: Let be an op- 

timal set of pure states. Then the decomposition 
{(^JV'i))}r=i =: representing the state p = 

S"=i is optimal. We now define a family of 

unitary (n -I- 1) x {n -\- l)-matrizes by U{t) := exp(tT), 
where t is a real parameter and T is the skew hermitian 
matrix defined by Tn+i,i = T^^n+i 
otherwise, i.e. 



Ci and Tij = 



T = 
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Selecting the first n columns of U{t) we get a right uni- 
tary (n-|- 1) X n-matrix U{t). Applying U{t) to the opti- 
mal decomposition according to eq. (||) we get 
a new decomposition {|0j(i))}j'=i of P which is defined 

by 



^2 " 



(5) 



1=1 



This decomposition can not have less entanglement than 
the decomposition which was assumed to be 

optimal, thus 

n n+1 

E^d^^)) <E^(i^j-w)) ^^^^ (6) 

holds true. At t=0, both sides of this inequality are equal. 
In the following, we will show that the first derivatives 
with respect to t of both sides are also equal at i = 0. 
Hence, the inequality must be fuUfiUed for the second 
derivatives at i = which will lead to eq (^. To proceed, 
we need the following lemma: 

Lemma: Let A{t) be a family of positive matrices, 
which depend differentiably on f € (a, b) . Then we have 



d , A{t) , 



Additionaly, if -^A{t) is zero at some point to and if A{t) 
is twice differentiable at to we have 



tv[A{to)" log2 



Ajto) 
tiAitoY 
(8) 



Note that eq. (0) is always well defined due to the pos- 
itivity and differentiability of A{t) whereas eq. (^) may 
not be well defined. 



Proof. Let A{t) = U {t)disig{X{t))W (t) with U{t) unitary. 



3 



Then we have 



d 



-( 5] A.(t)(log2 A. W - log2 5] A, W)) 



= -^AUi)(log2A.(t)-log2^A,(i)) 

i 

= -tr[^'(0(log2 
= - tr[A(t)' log2 



tr[^'(0(log2A(t)-log2 trA(t))] 
A{t) ^ 



tTA{t)'' 

In the last but one step, we have used that U'{t)W{t) + 
U{t)U'^'[t) = 0, that the trace is invariant under cychc 
permutations of the factors and that diag(A(t)) and 
log2[diag(A(i))] commute. To prove the second part of 
the lemma, we use the fact that if A' {to) is zero A -(to) is 
zero and we can always choose U{t) such that U'{tQ) — 0. 
This choice leads to 

^ -J2 A/'(io)(log2 A.(io) - log2 J2 ^j (^o)) 

i 3 

A{to) , 



-tr[A"(io)log2 



trA(to) 



having used U"{to)W{to) + U{to)W"{to) = 0. which 
holds true, because U'{to) = 0. □ 

Applying this lemma to Tri{t) := tr^ i — 

1, . . . , n, we get 

= -4(tr7f.(t)log2 ^'^'^ ' 



dt 



dt 

-tT^[{t) log2 



tTTTi{t)' 
n^{t) 



tTTTi{t) 



According to eq. (^, 7r^(i) is zero at i = 0. Therefore we 
can calculate the second derivative of £{\(j)i{t))) at t = 
using the second part of the lemma, which leads to 



dt^ 



t=0 



tr^",(0)log2 



MO) 
tr^,(0) 



Using eq. (|^) we can further evaluate the righthand side. 
Defining aij := tr^ we arrive at 



dt^ 



t=0 



\sV (J %- 



= 2 Re ( V c,Cj tr[ay log2 ) 

2 

n 



Re(^^CiCjtr[crjjlog2cr.ri]y (9) 



This expression is always well defined because 
tr Uij log2 Uii is well defined, even if ua has a non 
vanishing kernel. This can be easily seen by expressing 
I'f/'i) and in the Schmidt-basis of and then 

evaluating tr logj oa in this basis. 

To calculate the derivatives of £(|0„+i(i))) we use the 
fact, that E(c\f])) = |cp£(|?7)) which leads to 



-£(|0„+i(t)))| 



and 



1=1 





-o{e))\\ 



t=0 



dt 



(10) 



Note that the first and second derivative of 

^(Er=i c^lV-i) + C(i^))]|t^o ^^11 defined. 

As we have shown, the first derivatives of the right and 
left hand side of eq. m) are equal to 0. Therefore 



i=l 



(11) 



necessarily holds true, which in turn can be seen to be 
equivalent to eq. (^ by substituting eqs. (||) and ( |lO| ) 
into the right hand side as well as taking into account, 
that the left hand side of eq. ( pT| ) does not depend on t 
and therefore is zero. □ 

Our condition for the optimality of a set of pure states 
is a generalization of a result by Benatti and Narnhofer 

, who proved a similar condition for the special case of 
a decompositions consisting only of two vectors. 

It is easily seen, that it is sufficient to demand eq. (0) 
to be true for all Ci G C with ||c|| = 1 because we can 
cancel the square of the norm of the coefficient vector c 
on both sides of the inequality. Additionally, if the set 
{|V'i)}?=i is optimal, the sum on the right side of eq. (^) 
is already real. This is implied by the fact, that for two 
states lipk) and \ipi) from an optimal set, we have 

tr o-fci log2 o-fcfc = tr cr;fc log2 CTH , crfci := trs | V'fc)(V''i I ■ (12) 
To prove this equation we define the decompositions 

{i^.(^))}r=i with 

|0±(0)> = -^{coseii^k) ± sin 9\^i)) 
\if[B)) = ^(Tsin^lV^fc) +cos0|Vi)) and 



|0±(0)) = ^IV.) i + Kl 



4 



and {|?7.(0))}r^i with 



|ry±(0)) - ^(cos0|Vfc)±/sin0|V'i)) 
|77±(f?)> - ^(±/sin0|V'fe)+cos0|i/;,)) and 

For = those decompositions are identical to 
{^l^i)}"^]^ which is optimal by assumption. Hence 



> 



have to hold true, implying eq. (12) by using the lemma 
and some basic algebra. 

This theorem is usefuU to verify the optimality of a 
given decomposition or a given set of pure states, be- 
cause we only have to optimize a function over the unit 
sphere in C". In contrast, by using the straight forward 
condition 

n m n 

2—1 i — 1 2—1 

for the optimality of the decomposition {|'!/'i)}"=i which 
is easily derived from the definition of an optimal de- 
composition and eq. (|^), we had to optimize a function 
over the set of all right unitary m x n-matrices. In this 
case m must be at least [rankp]'^ with p = X]"=i 
which is the lowest known bound on the number of pure 
states in an optimal decomposition j^. Clearly this task 
is harder to perform. However, as there is nothing for 
free, the information we get in both cases is different: In 
doing the optimization over the right unitary matrizes 
we arrive at an optimal decomposition from which the 
entanglement of formation can be easily calculated. In 
doing the optimization over the unit sphere in C" , wc 
only get the information whether the decomposition was 
optimal or not. Therefore our condition is of no direct 
use in finding optimal decompositions and thus calculat- 
ing the entanglement of formation, but there is another 
important application, namely testing the additivity of 
the entanglement of formation. To perform this task, we 
have to check, whether the set of tensorproduct states 
{lip}) <8) IV'I)}"}^! is optimal, given two optimal sets of 
states {|V'i^)}iLi E^nd Numerically this task is 



easier to perform using our theorem, than using eq.(13) 
It might even be possible, that the new condition is help- 



Another information we get from our condition is an 
estimation of the entanglement of superpositions of states 
from an optimal set in terms of the entanglement of the 
states themselves and some overlap-terms. This can be 
seen by writing the diagonal terms in eq. (0) separately. 



i=l i=l 
n 

- Re ( ^ CiCj tr [cTy log2 an]) V Cj G C. 



Finally we conclude that in some sense locally optimal 
decompositions are already optimal: In deriving eq. (^) 
we have only used derivatives and the value itself at the 
point t = in eq. (|^) . Therefore it is sufficient to demand 
eq. to be true only in a small neighbourhood of t = 
0. This in turn is equivalent to demanding, that the 
decompositions {|'!/'i)}"=i h&s lower average entanglement 
than every decomposition, which we get from it by a 
right unitary (n + I) x n-matrix, which consists of the 
first n columns of a unitary (n -|- 1) x (n -|- l)-matrix 
in an arbitrary small neighbourhood of the identity, i.e. 
demanding local optimality. But since eq. (^) is also a 
sufficient condition for optimality, we conclude that local 
optimality implies optimality. 
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